Transmission properties of left-handed band-gap structures 
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We analyze transmission of electromagnetic waves through a periodic band-gap structure consist- 
ing of slabs of a left-handed metamaterial and air. Using the effective parameters of the metamaterial 
derived from its microscopic structure, we study, with the help of the transfer-matrix approach and 
by means of the finite-difference-time-domain numerical simulations, the transmission properties of 
such a left-handed photonic crystals in a novel type of band gap associated with the zero averaged 
refractive index. We demonstrate that the transmission can be made tunable by introducing defects, 
which allow to access selectively two different types of band gaps. 

PACS numbers: 42.70.Qs, 41.20.Jb, 78.20.-e 
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I. INTRODUCTION 

Materials with both negative dielectric permittivity 
and negative magnetic permeability were suggested the- 
oretically a long time ago 0, and they are termed left- 
handed materials because the wave vector and Poynt- 
ing vector lie in the opposite directions. Many unusual 
properties of such materials can be associated with their 
negative refractive index, as was demonstrated by several 
reliable experiments |^|^ and numerical finite-difference- 
time-domain (FDTD) simulations (see, e.g., Ref. Q). 

Multilayered structures that include materials with 
negative refraction can be considered as a sequence of 
the flat lenses that provide an optical cancellation of 
the layers with positive refractive index leading to ei- 
ther enhanced or suppressed transmission |5|, |6i, i3 • More 
importantly, a one-dimensional stack of layers with al- 
ternating positive and negative-refraction materials with 
zero averaged refractive index displays a novel type of 
the transmission band gap 01 S S Hll near the fre- 
quency where the condition <n> = is satisfied; such 
a novel band gap is quite different from a conventional 
Bragg reflection gap because it appears due to completely 
different physics of wave reflection. In particular, the 
periodic structures with zero averaged refractive index 
demonstrate a number of unique properties of the beam 
transmission observed as strong beam modification and 
reshaping being also insensitive to disorder that is 
symmetric in the random variable . 

In this paper, we study both transmission properties 
and defect-induced tunability of the left-handed photonic 
band gap structures created by alternating slabs of posi- 
tive and negative refractive index materials with an em- 
bedded defect taking into account realistic parameters, 
dispersion and losses of the left-handed material. We con- 
sider a band-gap structure schematically shown in Fig.^ 
First, we study the properties of the left-handed mate- 
rial as a composite structure made of arrays of wires and 
split-ring resonators (see the insert in Fig. 1) and de- 
rive the general results for the effective dielectric per- 
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FIG. 1: Schematic of a multilayered structure consisting of 
alternating metamaterial slabs and air. The inset shows the 
unit cell of the metamaterial structure. 



mittivity and magnetic permeability. Second, we study 
the transmission of electromagnetic waves through the 
layered structure consisting of alternating slabs of com- 
posite left-handed metamaterial using the calculated ef- 
fective parameters. We assume that the structure in- 
cludes a defect layer (see Fig. ^ that allows tunability 
of the transmission near the defect frequency. Using the 
transfer-matrix method, we describe the defect-induced 
localized states in such a structure and reveal that the de- 
fect modes may appear in different parameter regions and 
for both <n> ~ and Bragg scattering band gaps. De- 
pending on the defect parameters, the maximum trans- 
mission can be observed in all or just some spectral band 
gaps of the structure. We demonstrate that the frequency 
of the defect mode is less sensitive to manufacturing dis- 
order for the larger defect layer. Finally, we perform 
two-dimensional FDTD numerical simulations based on 
the microscopic parameters of the left-handed material 
and study the temporal evolution of the transmitted and 
reflected fields. 
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II. METAMATERIAL CHARACTERISTICS 

We assume that the left-handed metamaterial is cre- 
ated by a three-dimensional lattice of wires and split- 
ring resonators (SRRs), as shown in the inset of Fig. ^ 
Accordin g to the results of the derivation presented in 
Refs. mill, the main contribution to the effective di- 
electric permittivity of this structure is given by the 
wires, whereas the magnetic response is determined by 
SRRs. Although a three-dimensional lattice of wires 
contains closed circuits, we neglect their contribution to 
the magnetic permeability, because this effect is non- 
resonant, and, therefore, is weak. The effective dielectric 
permittivity can be obtained in the form [l^ 



Uj{uj - Z7e)' 



(1) 



where ujp ~ {c/d)[2TT/ In (d/rw)]^^^ is the effective plasma 
frequency, 7^ = /2aSlTi{d/rw), a is the wire conduc- 
tivity, S is the effective cross-section of a wire, S = irr^ , 
for 5 > r^,, and S ~ TrS(2r — S), for 5 < r^, where 
S = c/\/2n(TLj is the skin-layer thickness. 

To calculate the effective magnetic permeability of the 
lattice of SRRs we write its magnetization M in the form 
(see also Ref. [H) 



(2) 



where rim — i/d^ is the number of SRRs per unit cell, 
Rr is the SRR radius (see the insert in Fig.^, is the 
current in the SRR. We assume that SRR is an effec- 
tive oscillatory circuit with inductance L and resistance 
R of the wire, and capacity C of the SRR slot. In this 
circuit the electromotive force in this circuit due to an 
alternating magnetic field of the propagating wave. Un- 
der these assumptions, the evolution of the current in 
single SRR is governed by the equation 



L 



d'^Ir pdlr 1 ^ d£ 



(3) 



with 



£ = - 



7ri?2 dH' 
c dt 



where H' is the acting (microscopic) magnetic field, 
which differs from the average (macroscopic) magnetic 
field. We describe the SRR array as a system of mag- 
netic dipoles, which is valid when the number of SRRs in 
the volume A"^ is big enough, and use the Lorenz-Lorentz 
relation between the microscopic and macroscopic mag- 
netic fields [la 



H' = H+^M = B-^M. 
3 3 



(4) 



As a result, from Eqs. @~0 we obtain the effective 
magnetic permeability of the structure in the form. 



(w) = 1 + 



Lu^ -u^ [1 + F/3) + iuj"f' 



(5) 
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FIG. 2: Frequency dependence of the real part of the dielec- 
tric permittivity Re(e) (solid), and the real part of magnetic 
permeability Re(/i) (dashed). 



where F = 2TTnm{Tr Rl/ cf / L , ujI = 1/LC, and 7 = R/L. 
Inductance L, resistance R, and capacitance C are given 
by the following results (see, e.g., Ref. |0|), 



L = 



In 



R = 



2TlRr 
O Sr 



c = 



4:TTdo 



where r is the radius of the wire that makes a ring, a 
is conductivity of the wire, Sr is the effective area of 
the cross-section of SRR wire defined similar to that of 
a straight wire, and dg is the size of the SRR slot. We 
note, that the result for C should hold provided dg <C r. 

Taking the parameters of a metallic composite as d = 1 
cm, r„ ~ 0.05 cm, Rr — 0.2 cm, r = 0.02 cm, dg = 
10"'^ cm, and its conductivity as ti = 2 • 10^^ s~^, we 
calculate the effective frequency dependencies of e(w) and 
fi{uj) according to Eqs. Q and lO, respectively, and show 
these dependencies in Fig. |5] The resonance frequency 
appears near 5.82 GHz, and the region of simultaneously 
negative e and fj, is between 5.82 GHz and 5.96 GHz. 
The imaginary part of the magnetic permeability, which 
determines effective losses in a left-handed material, is 
larger near the resonance. 



III. TRANSMISSION AND DEFECT MODES 

We consider now a band gap structure created by seven 
alternating left-handed and dielectric slab pairs, as shown 
in Fig.^ The number of slabs is chosen to keep losses in 
the structure at a reasonably low level, still having visible 
the effects of periodicity. We assume that the periodic 
structure is created by slabs of the width a, made of the 
left-handed composite described above, and separated by 
air (eq — 1, fia = !)■ The middle layer of the left-handed 
material is assumed to have a different thickness, b = 
a(l -I- A), in order to describe a structural defect. 

To study the transmission characteristics of such 
a band-gap structure, we consider a normal-incidence 
plane wave for the scalar electric field described by the 



Helmholtz-type equation, 



dz2 



^e(z)^(z) - 



dfi d 



^{z) dz dz 



(6) 



where e{z) and fi{z) are the dielectric permittivity and 
magnetic permeabihty in the whole space. 

Before analyzing the transmission properties of a finite 
layered structure, first we study the corresponding infi- 
nite structure without defects and calculate its band gap 
spectrum. In an infinite periodic structure, propagat- 
ing waves in the form of the Floquet-Bloch modes satisfy 
the periodicity condition, E{z + 2a) = E{z) exp{2i a Kh), 
where Ki, is the Bloch wavenumber. The values of Kh are 
found as solutions of the standard eigenvalue equation for 
a two-layered periodic structure (see, e.g., Ref. fTHl. 



2 cos(i^f,2a) — 2 cos [{kr + ki)a] — 

Pr_ , Pi 
Pi 



Pr 



+ 2 ) • sin {kro) ■ sin (fcja). 



(7) 



where Pr — I, pi — y'e/ p,, kr = lo/c and ki = lj / c^jeji are 
the wavevectors in air and left-handed slabs, respectively. 
For real values of Ki,, the Bloch waves are propagating; 
complex values of Ki, indicate the presence of band gaps, 
where the wave propagation is suppressed. The spectral 
gaps appear when the argument of the cosine function 
in Eq. Q becomes the whole multiple of tt, and no real 
solutions for Kh exist. These gaps are usually termed 
as Bragg gaps. The presence of the left-handed material 
in the structure makes it possible for the argument to 
vanish, so that the wave propagation becomes prohibited 
in this case as well, thus creating an additional band gap, 
which do not exist in conventional periodic structures. As 
a matter of fact, the condition \kr\ = \ki\ corresponds to 
the zero average refractive index, <n> = 0, as discussed 
in Refs. 0, However, the inherent feature 

of the left-handed materials is their frequency dispersion, 
so that the condition |fc,.| = \ki\ defines a characteristic 
frequency w* at which the indices of refraction in both 
the media compensate each other. In a sharp contrast 
to the conventional Bragg reflection gaps, the position of 
this additional <n> = gap in the frequency domain 
does not depend on the optical period of the structure. 

For the parameters of the left-handed media described 
above, the frequency lo* at which the average refractive 
index vanishes is found as w* « 27r • 5.8636 x 10^ s~^. 
Importantly, the transmission coefficient calculated at 
uj — lli* for the seven-layer structure shows a characteris- 
tic resonant dependence as a function of the normalized 
slab thickness a/A, where A = 2t:c/uj*, as shown in Fig.|31 
The transmission maxima appear in the <n> = gap, 
when the slab thickness a coincides with a whole multiple 
of a half of the wavelength. The width of the transmission 
peaks decreases with the increase of the number of lay- 
ers in the structure. The transmission maxima decrease 
with increasing thickness of the structure due to losses in 
the left-handed material which become larger for thicker 
slabs. One of the interesting features of the <n> = 




FIG. 3: Transmission coefficient of a finite periodic structure 
with seven layers of left-handed material vs. the normalized 
slab thickness a/A, where A = 2nc/Lj* and the frequency uu = 
uj* corresponds to the condition <n> = 0. 
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FIG. 4: Transmission coefficient of the left-handed band-gap 
structure vs. the wave frequency, (a) The structure with the 
period a = 0.25A, without (solid) and with (dashed) defect 
layer (A = —0.8). (b) The structure with the period a = 
1.25A without (solid) and with (dashed) defect layer (A = 
-0.6). 



gap is that the transmission coefficient can vanish even 
for very small values of the slab thickness. This prop- 
erty can be employed to create elective mirrors in the 
microwave region operating in this novel <n> = gap 
which can be effectively thinner than the wavelength of 
electromagnetic waves. 

The transmission coefficient is shown in Fig. 01a,b) as 
a function of the frequency, for two structures with dif- 
ferent values of the period a. For the quarter-wavelength 
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FIG. 5: Frequency spectrum of the defect modes as a function 
of tiie normalized defect size A in tlie left-lianded band-gap 
structure with the period a = 1.25A. 
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FIG. 6: Numerical FDTD simulation results showing relax- 
ation processes in a band-gap structure with a defect. Solid - 
incident energy flow, dashed - transmitted energy flow, dotted 
- reflected energy flow. Parameters are: a = 0.25A, A = —0.8. 
(a) Defect mode is not excited, w = 27r5.86 x 10^ s~^; (b) De- 
fect mode is excited, u) = 27r5.878 x 10® . 



slabs [see Fig. EI^)], the only visible band gap is the 
<n> = gap centered near lo* where both e and [i are 
negative. When the structure has a defect, the trans- 
mission peak associated with the excitation of the de- 
fect mode appears inside the <n> = gap as shown by 
the dashed curve. For the structure with thicker slabs, 
e.g. for a = 1.25A [see Fig. ^b)] the <n> = gap 
becomes narrower but it remains centered near the fre- 
quency uj* . The transmission coefficient of this second 
band-gap structure shows, in addition to the <n> ~ 
gap, two Bragg scattering gaps. Due to the increased 
losses in this second band-gap structure, where slabs are 
thicker than in the structure corresponding to Fig. 2Ja), 
the effects of the resonant transmission at the defect 
mode become less visible. Moreover, for the parameters 
we consider here the defect mode appears only in the 
<n> — gap, whereas it does not appear in the Bragg 
gaps. For larger slab thickness, higher-order Bragg gaps 
may appear in the frequency range where the composite 
material is left-handed. 

In Fig.|Sl we show the frequency spectrum of the defect 
modes for the structure with a = 1.25A as a function of 
the normalized defect size A. We notice a number of 
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FIG. 7: Numerical FDTD simulation results for the pulse 
scattering by a periodic structure with defect. Solid - inci- 
dent energy flow, dashed - transmitted energy flow, dotted - 
reflected energy flow. Parameters are: a = 0.25A, A = —0.8. 
(a) Defect mode is not excited, u) — 27r5.86 x ICp s~^; (b) 
Defect mode is excited, u — 27r5.878 x 10® . 



important features: (i) the defect modes do not always 
appear simultaneously in all gaps, and (it) the slope of the 
curves in Fig. [31 decreases with the growth of the defect 
thickness. As a result, the eigenfrequencies of the modes 
introduced by thicker defect layers can be more stable 
to disorder introduced by manufacturing. These novel 
features seem to be important for tunable properties of 
the layered structures because the different modes allow 
to access different types of band gaps. 



IV. NUMERICAL FDTD SIMULATIONS 

In order to analyze the temporal evolution of the trans- 
mitted fields and the beam scattering under realistic con- 
ditions, we perform two-dimensional FDTD numerical 
simulations of the beam propagation through the left- 
handed periodic structure of seven layers with a defect. 
We consider TM-polarized Gaussian beam of the width 
20A propagating towards the structure with the period 
a = 0.25A; such a structure corresponds to the trans- 
mission coefficient shown in Fig. 01 a) by a dashed line. 
First, we choose the frequency of the incident field in 
the <n> = gap. The temporal evolution of the en- 
ergy flows (integrated over the transverse dimension) for 
the incident, transmitted, and reflected waves is shown 
in Fig. ini^a), clearly indicating that the transmission 
through such a structure is negligible. When the fre- 
quency of the incident field coincides with that of the de- 
fect mode, a significant amount of the energy is transmit- 
ted through the structure [see Fig.lBfb)]. The relaxation 
time of the beam transmission through the structure is 
estimated as 10'^ periods (approximately 170 ns). 

Results of FDTD simulations for the pulse scattering 
from the structure with a — 0.25 A are shown in Figs. [7| 
(a,b) as the temporal dependence of the incident, re- 
flected and transmitted energy flows. One can clearly 
see significant transmission, when the carrier frequency 
of the pulse coincides with the eigen frequency of the 
defect mode [see Figs. [T] (a,b)]. 
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FIG. 8: Results of the numerical FDTD simulations for the 
amplitude of the magnetic field in a two-dimensional structure 
(natural logarithm scale). Boxes show positions of the left- 
handed slabs, a = 0.25A, A = —0.8. (a) Defect mode is not 
excited, lo = 27r5.86 x 10^ s~^; (b) Defect mode is excited, 
uj = 27r5.878 x 10^ 



An example of the field distribution in the structure 
with the slab size a — 0.25A is shown in Figs. |S1 (a,b) 
for two different regimes. In Fig. Eta), the frequency 
corresponds to low transmission in Fig. ^a), when no 
defect mode is excited. Figure |HIb) demonstrates the 
field distribution in the structure with an excited defect 
mode and enhanced transmission. 



V. CONCLUSIONS 

We have studied the transmission properties of peri- 
odic structures made of a left-handed metamaterial and 
air. Using realistic parameters of the metamaterial de- 
rived from the microscopic approach, we have calculated 
the band-gap spectrum of an infinite one-dimensional 
structure with and without defects, demonstrating the 
existence of band gaps of two different types, the conven- 
tional Bragg scattering gaps and a novel <n> = gap. 
We have analyzed the properties of the defect modes in 
a finite periodic structure with a structural defect and 
demonstrated that, depending on the defect size, the de- 
fect modes appear in all or just some of the band gaps al- 
lowing to access different gaps selectively. In addition, we 
have performed two-dimensional numerical FDTD simu- 
lations of the propagation of electromagnetic waves in 
such structures and have studied the temporal dynamics 
of the beam transmission and reflection. We have demon- 
strated that the excitation of defect modes can enhance 
substantially the wave transmission through the struc- 
ture. 

The authors thank Michael Feise for collaboration, and 
acknowledge a partial support from the Australian Re- 
search Council. 
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